o 

(N 



l> 



> 

o 



Existence and uniqueness of very singular 

solutions for a fast diffusion equation 

with gradient absorption* 



Razvan Gabriel lagar ]^ 
Philippe Laurengot ^ 



Abstract 



Existence and uniqueness of radially symmetric self-similar very singular solutions 
P^ ■ are proved for the singular diffusion equation with gradient absorption 
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1 Introduction 

The singular diffusion equation with gradient absorption 

dtu- ApU+\Vu\'^ = 0, (i,x) gQoo := (0,00) xM^, (1.1) 

with p € (1,2) and q > features a singular diffusion term and an absorption term 
depending solely on the gradient. Taken apart, these two terms lead to two completely 
different behaviors for large times so that, putting them together, a competition between 
them is expected. In fact, we study recently in [10] qualitative properties and decay 
estimates for nonnegative and bounded solutions to (jl.ip and identify ranges of exponents 
p S (1,2) and q > with different behaviors. In particular, given p E [2N/[N + 1),2), 
if g > Q* := p — {N/[N + 1)), the diffusion term dominates for large times while it is 
the absorption term that dominates when q ^ (0,p — 1) leading to finite time extinction. 
Finite time extinction also occurs for q (z {p — l,p/2) but is expected to be of a different 
nature, some influence of the diffusion term persisting near the extinction time. Finally, 
when q G (p/2,g*), the solutions to (jl.ip with initial data decaying sufficiently rapidly 
at infinity decay to zero at a faster algebraic rate than the one that would result from 
the diffusion alone, a feature which reveals an interplay between diffusion and absorption. 
By analogy with the existing literature on related problems (see, e.g., [21 |8l [23] and the 
references therein), the large time behavior in that case is expected to be described by a 
particular self-similar solution to (jl.ip which is called a very singular solution. Recall that 
a very singular solution to a partial differential equation is a solution U in Qoo (in a weak 
or classical sense) such that 

lim sup{C/(t,3;)} = 0, (1.2) 



*^0|x|>e 



and 



lim / U{t,x) dx = 00, (1-3) 



*^0-'|x|<e 



for any e > [5]- In other words, the initial condition for U is zero in M \ {0} and 
it has a stronger singularity at the origin x = than the Dirac mass. Recall that a 
solution to a partial differential equation is usually referred to as a singular solution if 
it satisfies (|1.2p . the main examples being the so-called fundamental solutions, that is, 
solutions having a Dirac mass as initial condition, and the very singular solutions. The 
existence and non-existence of singular solutions have been thoroughly studied for diffusion 
equations with a zero-order absorption term and we refer to, e.g., [5l [12l [HI [151 [13 [E] 
for dtu - Au"" + u"? = 0, m > 0, and m [13 [IH] for dtu - ApU + ui = 0, p > 1. The case 
of diffusion equations with an absorption term depending solely on the gradient has been 
investigated more recently, see, e.g., [Tl[3l[20] for dtu — Au™ + |Vti|'^ = 0, m > 1, and [21] 
for dtu - ApU + |Vn|'? = 0, p > 2. 

In this paper we focus on the singular diffusion equation with gradient absorption (jl.ip 
for the particular range of exponents p and q for which very singular solutions are likely 
to exist, namely: 

2N p N , ,, 

as already mentioned. Owing to the homogeneity of (jl.ip . we actually look for a (forward) 
self-similar and radially symmetric very singular solution u to p.ip of the form 

n(t,x)=t-"/(|x|t-^), (t,x)GQoo, 



for some profile / and exponents a and /3 to be determined. Inserting this ansatz in (jl.ip 
gives the values of a and /3 

P-Q a g-P+1 /. KN 

2g - p 2g - p 

and implies that the profile / = /(r), r = \x\t~^, is a solution of the ordinary differential 
equation 

(l/'r'/O'W + ^^(l/'r'/')(0 + «/(r) + /3r/'(r) - |/'(r)r = , r > 0, (1.6) 

with /'(O) = 0. The previous conditions (|1.'2|) and (|1.3|) become 

lim r(P-5)/('?-P+^)/W = 0, limt^^~" / /(r)r^-i dr = cx), (1.7) 

the second one being satisfied if / G -^^(0, oo; r^~^dr) and a — Nf3 > 0, that is g < g*. 
We then prove the following result: 



Theorem 1.1. Assume that p and q satisfy (|1.4p . There is a unique nonnegative solution 
f to ()1.6p which satisfies /'(O) = and 

lim ^(p-^VC^-P+i) j(^) = 0. (1.8) 

r— i-oo 

In fact, there is a positive constant w* defined in ()2.19p below such that 

lim rP/^^'P'^ f{r) = w* . 

r— >oo 

Let us first mention that Theorem 1.1 is given in [22j under the additional restriction that 
g > 1 besides the constraints (|1.4p on p and q. However, the proof given there does not 
seem to apply to the case q G (p/2, 1] and some new ideas have to be introduced which 
actually work for the whole range defined in (|1.4p . The proof given below is thus done for 
p and q satisfying ()1.4p . 

As usual, Theorem 11.11 is a consequence of a detailed study of the initial value prob- 
lem associated to (|1.6p . that is, we consider the solution /(.;a) to (jl.6p with initial data 
/(0;a) = a > and /'(0;a) = 0. We first establish the well-posedness of this problem 
at the beginning of Section [2] together with some basic properties of its solutions in Sec- 
tion l2.1l In particular, we show that, if a > is such that /(.; a) is positive in (0, cxd), then 
!/'(.; a) I is controlled by /(.; a)^'^ as a consequence of sharp gradient estimates established 
in |10j . This turns out to be a cornerstone of the proof as it allows us to transfer some 
asymptotic properties of /(.; o) to /'(.; o). We next prove in Section [2.21 that a monotonic- 
ity property with respect to the initial condition a is enjoyed by the solutions /(.; a) under 
suitable conditions. Such a monotonicity property is also true for the singular diffusion 
equation with zero-order absorption studied in [6] . We then split the range of initial con- 
ditions a G (0, oo) into three disjoint sets according to the behavior of the derivative of 
r I—)- rPi^'^~P> f{r;a) and characterize these three sets in Sections 12.31 12.41 and 12.51 This 
study then guarantees the existence of at least one profile / satisfying the properties listed 
in Theorem 11.11 We finally borrow an argument from [6] to prove the uniqueness of this 
solution. 



2 Self-similar and radially symmetric very singular solutions 

In order to construct a solution to p.6|) satisfying (|1.8p . we use a shooting method which 
leads us to the following initial value problem: 

(l/TVO'W + ^^^(l/'r'/OW + «/« + Prfir) - |/'(r)r = 0, 

(2.1) 

I /(O) = a, /'(O) = 0, 

where a > is an arbitrary positive real number and the condition /'(O) = follows from 
the radial symmetry and expected smoothness of very singular solutions. Introducing 
F := — 1/'|^~^/', we observe that (|2.ip also reads 

/'(r) = -|F(r)|(2-p)/(p-i)F(r), 

F'(r) + ^l^F(r) = af(r) - /3r|F(r)|(2-P)/(P-i)F(r) - \F(rW/'-P~^\ (2.2) 

r 

I /(O) = a, g(0) = 0. 

Since q > p — 1 and (2 — p)/{p — 1) > 0, the right-hand side of (j2.2|) is locally Lipschitz 
continuous. The term involving {N — l)/r being handled as usual, there is a unique 
maximal C^-smooth solution (f{.;a),F(.;a)) to (j2.2p . Owing to the continuity of /(.;a) 
and the positivity of a, f{.;a) is clearly positive in a right-neighborhood of r = and we 
define 

i?(a) :=inf{r > : /(r;a) = 0} > 0. (2.3) 

In the sequel, where there is no risk of confusion, we will omit a from the notation and let 

/ = /(•;«)• 

2.1 Basic properties of /(.;a) 

We first prove some basic properties of the profile / = /(.; a). 

Lemma 2.1. Let a > 0. We have f'{r;a) < for any r € (0, i?(a)) and \f'{r;a)\ < (aa)^^'^ 
in (0, i?(a)). Moreover, if R{a) = oo, then lim f{r;a) = lim f'{r;a) = 0. 

r— >oo r— >oo 

Proof. It readily follows from ()2.2p that F'{0) = aa/N > so that there is 5 > such that 
/'(r) < for r € (0,5). Introducing tq := inf {r G (0, i?(a)) : /'(r) = 0}, we assume for 
contradiction that tq < R{a). Then, on the one hand, F{rQ) = /'(ro) = and we deduce 
from ([221) that F'(ro) = a/(ro) > 0. On the other hand, F{r) > = F{ro) for r G (0,ro) 
which implies that -F'(ro) < 0, whence a contradiction. Consequently, rg > i?(a) and 
/' < in (0, R{a)). Consider now R G (0, R{a)) and let rm be a point of minimum of /' in 
[0, R]. Clearly, r^ > and either r^ G (0, R) and f"{rm) = or r^ = R and f"{rm) < 0. 
In both cases, it follows from (j2.ip that \f'{rm)\'^ — af{rm) < 0, whence 

If'irW < l/'(^m)r < «/(^m) < «/(0) = aa, re[0,R]. 

Assume now that R{a) = oo. Since / is decreasing, there exists / := lim /(r) > 0. 

r— >oo 

Defining the following energy: 

E{r):=P^\f'{r)\P + ^f{r)^ r > 0, 
p 2 



it follows from (12.111 that 



E'{r) = {p-l){\f'r'f'n{r)+a{ff'){r) 

= -^^^^l/'(r)r-/3r|/'(r)|2-|/'(r)r+^<0. (2.4) 

r 

Hence E is decreasing and positive in (0, cxj), and has thus a limit as r — > oo. This property 
ensures that /' has a limit as r — > oo while (j2.4p implies that /' € ^'"^^{0, oo). Combining 
these two facts entails that f'{r) — > as r ^> oo and so does F. It then follows from ()2.3p 
that F'{r) — > al as r —^ oo. Assume for contradiction that / > 0. Then F'{r) > al/2 for 
r large enough whence F{r) > air /A for r large enough. Therefore F{r) — > oo as r ^- oo 
and a contradiction. We have thus shown that / = 0. □ 

Another useful result is the following expansion of /(.; a) near the origin: 

Lemma 2.2. Let a > 0. There exist three positive constants Ci, C2, and C3 such that, as 
r — ^ 0, 






and 



jI (r^-V'(OrV'(r)) = -aa+ (^)'^^'"^''/(-i) +o(W(-i)). 



r^ 1 (;;^^ 

Integrating once, we obtain the following expansion for /': 



;^(,) = _(-)^/(-^l/(.-l) 



^ q + N{p 



1 f^V'^"^+'^/^^-'^(.+i)/(p-i) +of^{.+i)/b-i)) 

(p-l)ViVy' ^"^ >■ 



(2.5) 



(2.6) 



•'^ ' ^ (jo-l)a ViV/ 

(2 + g-p)C2 /aa\ (^-P+2)/(P-i) (p+,)/(p_i) . (p+,)/(p_i). 
(p-l)a VAf/ ^ ^ ^' 

Proof. We start with the fact that (|/T~^/')'(0) = -(aa/A^), so that, as r ^ 0, 

After integration, we obtain 

/(r) = a - ^^ (^)'^^'"'^ -^/^^"'^ + o(r^/(^-^)). (2-8) 

Since (j2.ip also reads 

d_ 
dr 

we infer from ()2.7p and ()2.8p that, as r ^- 
1 d 



(r^-i|/'(r)r-V'(r)) = r^"^ (|/'(r)r - /^r/'(r) - a/(r)) , (2.9) 



(2.10) 



Another integration gives that, as r ^- 0, 

^ ^^^ (211) 

Inserting ()2.10p and ()2.1ip in ()2.9p then leads to ()2.5p after lengthy but simple computa- 
tions, with the constants 

^ P-1 ^ P-1 ^ _ (P-I)g 

^1 — ; ^2 — 7 ; 77 -^TTT TTT' ^3 



P ' (p + g)(g + iV(p-l))' 2p2(p + Ar(p_l))(2g-p)- 

Concerning daf, we first observe that differentiating (j2.9p with respect to a gives 

-/3r^9,/'(r) - ar^-iaa/(r) . (2.12) 
Setting Y'{r) := q\f'{r)\'i-'P f'{r)/{p - 1) for r G [0,i?(o)) with y(0) = 0, we have 

i. (r^- V'(r)r^9./'(r)e-^«) = - [Pr^daf{r) + ar^-^a,/(r)] ^ . 

We next use the properties daf{0;a) = 1, daf'{0;a) = 0, and Y'{r) — > as r — >■ (the 
latter being true since q > p — 1) to conclude that 



as r ^ 0. Integrating this identity, we infer from (|2.5p that 

■^ ^ ^ Af(p-i) viv; ^ ^ 

Integrating once more with respect to r gives an expansion of daf{r) as r — > 0. Insert- 
ing these expansions in (|2.12p and arguing as in the proof of (|2.5p give (|2.6p after some 
computations. □ 

We now prove that, if a > is such that R{a) = oo, then |/'(.;a)| is controlled by 
/(.;a)2'P. To this end, we first check that /(.;a) can be associated to a solution to (jl.ip 
which turns out to be a viscosity solution in the sense of Definition lA.ll below. This then 
allows us to apply the optimal gradient estimates obtained in |10j . 



Lemma 2.3. Let a > be such that R(a) = oo and e > 0. Setting 

Ue{t,x) := {t + eyf{\x\{t + e)-'^;a), {t,x) G [0, oo) x M^ , 

the function Us is a viscosity solution to (jl.ip with initial condition C4(0) defined by 
UeiO,x) ■.= e-'^f{\x\e-^;a) forxeR^. 



Proof. Let us first observe that, owing to ()2.7p - ()2.9p . we have as r — t- 

{p-l)f"{r) = - (^)'^^^"'V(2-P)/(P-1) +o(r(2-rt/(P-i)). (2.13) 

Since p < 2 and / clearly belongs to C^((0, oo)), we deduce from ()2.13p that / G C^([0, oo)) 
with /"(O) = 0. Consequently, Ue belongs to C'^{[0,qo) x R''^) and is clearly a classical 
solution to (jl.ip in (0, oo) x M \ {0} where VC/g does not vanish. Then Us is a viscosity 
solution to (jl.ip in (0, oo) x R^ \{0} as it obviously satisfies Definition lA.il for all (to, a^o) € 
(0,oo) xM^\{0}. 

Consider now tQ > and xq = 0. We first prove that Us is a viscosity subsolution. To 
this end, let V € ^ be such that Us{to,0) = ip{tQ,0) and Us{t,x) < '4^{t,x) for any (t,x) € 
Qoo\{(iO) 0)}, the set A of admissible comparison functions being defined in the Appendix. 
Since U^- and ip are both C^'^-smooth, this property implies that Wtpito, 0) = VC/e(io, 0) = 
and dtip{to,0) = dtUs{to,0) = —aa{to + e)"+^ < 0, so that Definition lA.il is satisfied. 

We next prove that Us is a viscosity supersolution, that is, — C/g is a viscosity subsolution. 
Let ip G A he such that —Us{to,0) = ijj{tQ,0) and —Us{t,x) < ip{t,x) for any (t,x) € 
Qoo \ {(^0)0)}. Since Us and ip are both C"^-smooth, we have Vip{tQ,0) = VUs{to,0) = 
and 9t^(io,0) = -dtUs{to,0). Therefore, for x € M^, 

tp{to,x) - i^{to,0) > -Us{to,x) + Us{to,0) = -{to + e)-'' [f{\x\ito + ey^;a) - f{0;a) 
and it follows from (|2.5p that, as x — >■ 0, 

V^(to,x) - V(io,0) > Ci(to + e)-" (^y^^^"'^ (j^|(to + e)-/3)^/^^"'Vo(|xr/(P-i)). 

However, since ip £ A, there is ^ G H such that, as x — > 0, ^(|x|) > ip{to,x) — ^(to,0). 
Combining the above two inequalities leads us to 

liminf y. \, >0, 

which contradicts ()A.2p . This situation thus cannot occur and the proof is complete. □ 

As a consequence of Lemma 12.31 and since p > pc, we may apply the gradient estimates 
proved for solutions to (|l.ip in |1U| Theorems 1.2 and 1.3]. 

Proposition 2.4. Let a > be such that R(a) = oo. Then f{.;a) satisfies 

|/'(r;a)|<C4/(r;a)2/P, r>0, (2.14) 

for some constant C4, depending only on N , p, q and a. 

Proof. Fix e > 0. According to Lemma 12.31 the function Us defined by Us{t,x) = 
{t + e)-"/(|x|(t + e)"^;a) for (i,x) G [0,oo) x M^ is a viscosity solution to ([II]) with 
a nonnegative initial condition which belongs to Ty^'°°(M^) by Lemma |2. 11 
If q > 1, we infer from [10| Theorem 1.2] that there is a positive constant C depending 
only on A^ and p such that 

VC/,-(2-p)/P(i^ x) < C t-^/P , (t, x) G Qoo. 



In terms of /(.; a), we obtain 



n\x\{t+er^) 



<C 



t + e 



i/p 



f{\x\{t + e)-^f/P, {t,x)eQ. 



Letting e — > and choosing t = 1 give (j2.14p . 

If g G (p/2,1), we infer from [101 Theorem 1.3] that there is a positive constant C 
depending only on A^, p, and q such that 

In terms of /(.; a), we obtain 

/'(|x|(t + e)-'^)| <C {t + ey/P (ai/"P(t + 2e)-i/P + t-^/p) /(|x|(t + e)-^)^/^ 
for {t,x) e Qoo- Letting e ^- and choosing t = 1 give (|2.14p . □ 

2.2 Decay rates and monotonicity 

We now study the possible decay rates as r — > cx) of the profile /(.; a) (when R{a) = oo). 
Since we expect an algebraic decay, we make the following ansatz: 



/(r; a) ~ Cr'-\ /'(r; a) ~ 0(1 - 7)r"^ 



as r 



oo, 



for some 7 > 1. Inserting this ansatz in (j2.1|) . we easily find that there are only two 
possibilities: 7 = 2/(2 —p) and 7 = l/{q — p+ 1). In the former case, /(r; a) ~ Cr~^'^'^~^> 
as r — 7- 00, while f{r;a) ~ C'r"''^"'^^' '''?~'P+^) as r ^ 00 in the latter. Observing that 
{v ~ q)/{q ~ P + 1) < p/(2 — p) as g > p/2, the solutions to ()2.ip decaying with the first 
rate are called fast orbits while those decaying with the second one are called slow orbits 
[5]. Complying with ()1.2p requires / to be a fast orbit, and we now proceed to show the 
existence of such solutions to (j2.ip . 
To this end, following [6l[22], we introduce the new unknown function 

w{r) = w{r; a) := rP/^'^-P^ f{r; a), (r, a) G [0, R{a)) x (0, 00), 

which is a solution of the following differential equation: 

{p - l)r'^w"{r) + [N -I- 2/i(p - l)]rw'{r) + /x(/i - N)w{r) 

_l_ |j-y;'(7-) _ ^w{r)\ ~P [(a — (i^)w{r) + /3rw'{r) — r'^\rw'{r) — fiw{r)\^~^ 

for r G [0,R{a)), where 







(2.15) 



^l:-- 



p 



> N, rj:-- 



2q — p 



< 0, a — /3fi 



<0. 



p 'I — p 2 — p 

We end the preliminary results with the following monotonicity property with respect to 
the initial value a which will be very useful in the uniqueness proof later on. We define a 
linear differential operator La by: 

La{^){r) := {p - l)rV"(r) + [A^ - 1 - 2fi{p - l)]r^'{r) + fi{fi - N^r) 

+ (2 - p){\W\-PW){r){rip'{r) - fiip{r)) [{a - Pfi)w{r) + /3™'(r) - r''|VF(r)|«] 

+ \W{r)\'^-P [{a - PuMr) + I3rip\r) - qr''{\W\'^-^W)ir){r^p'{r) - Mr))] , 

(2.16) 



where W{r) := rw'{r)—n'w{r) and w{r) = w{r; a) for r € [0, R{a)) in ()2.16p . Before stating 
the monotonicity lemma, we gather some properties of La in the fohowing two results. 

Lemma 2.5. Given a > 0, we have 

La{daUj{.;a)) =0 and La{rw\.;a)) < in {0,R{a)). (2.17) 

Moreover, for I € {0, fia], there exists a small interval {0,si{a)) such that 

fiawa{r;a) > rw'{r;a) for r e {0,s^a{a)) , 
lwa{r;a) < rw'{r;a) for I G {0, ij,a) and r€(0, s;(a)). 

Proof. To simplify notations, we set Wa := daw{.;a) and fa '■= daf{-]a). Differentiating 
the ordinary differential equation ()2.15p with respect to the parameter a, we readily obtain 
that La{wa) = in (0, R{a)). Next, differentiating (j2.15p with respect to r and multiplying 
the resulting identity by r give, by a straightforward computation 

La{rw') = rjr'^lrw'ir) - /xu;(r)|'?~P+2 < g in (0,i?(a)). 

We now prove the last two assertions and consider / € (0, /.m]. Since Wa{r) = r^fa{r) and 
rw'{r) = r^^{rf'{r) + iJ.f{r)), we have lwa{r) — rw'{r) = r>^ (Ifai^) — rf'{r) — nf{r)) and 
use Lemma 12.21 to identify the behavior of this function in a right-neighborhood of r = 0. 
Indeed, we infer from Lemma [2^21 (|2.10p . and the definitions of Ci and C2 that, as r — > 0, 



lfa{r) - rf'{r) - fif{r) 

p/(p-i) 



^ ,'aa\i/(p-i) 
l-Ha-Ci( — 



I p 



^ /aa\('?-P+2)/(p-i) 
+ o(r^'^~^'^^^^^~^h . 



{p — l)a p — 1 

{2 + q-p)l (p + q) 
{p — l)a p — \ 



J.I- 



/^ 



j-(p+'?)/(p-i) 



On the one hand, if / < ^a, we have lfa{r) — rf'[r) — fif{r) ~ / — /ua < as r —t- 0, which 
implies that this function is negative in (0, s/(a)) for some sufficiently small si{a) > 0. On 
the other hand, if I = fia, we have 

I p fi p 

l-fia = 0, /x = — - fi = 

[p — l)a p—1 [p — 1) p— 1 



and 



{2 + q-p)l {p + q) fi{2 + q-p) p + q {2q - p){p - I) 

fi = — H = > 0, 



(p — l)a p—1 p—1 p—1 2—p 

hence fiafa{r) — rf'{r) — fif{r) > for r > sufficiently small. □ 

Lemma 2.6 (Comparison principle). Let a > 0, ri S (0, -R(a)), and r2 G (ri,i?(a)) 
and assume that w'{.;a) > in [ri,r2]. Then, for any function h G C^(ri,r2) satisfying 
h{ri) = h{r2) = and La{h) >0 in (ri,r2), we have h < in (ri,r2). 

Proof. Lemma [2.61 is a variant of the classical maximum principle [HI Corollary 3.2] and is 
proved in [U p. 48]. The proof relies on the existence of a function 93 > on [ri,r2] such 
that La{^) < 0. Here this is satisfied by (p{r) = rw'(r;a) in view of Lemma [231 and the 
assumption on !(;'(.; a). [] 



We are now in a position to state and prove the monotonicity result with respect to a. 

Lemma 2.7 (Monotonicity lemma). Let a > and suppose that w'{.;a) > in a finite 
interval (0, rg) with tq < R{a). Then we have fiadaw{r;a) > rw'{r;a) for r G (0,ro) and 
daw{.;a) > on [0,ro]. 

Proof. We adapt the proof of [Q Lemma 2.2]. To simplify notations, we again set Wa '■= 
davu{.;a) and fa := daf{-;a) and first show that nawa{r) > rw'{r) for r E (0,ro). Assume 
for contradiction that there exists some point r2 € (0, tq) such that fiaWa{r) — rw'{r) > 
for r € (0, r2) and fiawa{r2) — r2w'{r2) = at r = r2. Then, given / G (0,/ia), we have 
lwa{r2) — r2w'{r2) < and, since lwa{r) — rw'{r) < for r G (0,Si(a)) by Lemma \2.h\ 
there is I G (0, //a) sufficiently close to ^a such that we have lwa{ri) —riw'{ri) = lwa{rz) — 
r2,w'{rj,) = and lwa{r) — rw'{r) > for r G {ri^r-^) for some < ri < ra < r2. We 
are then in a position to apply Lemma 12.61 to the function h{r) := lwa{r) — rw'[r) in 
the interval [ri,r3] C (0,r2), recalling that La{h) = lLa{wa) — La{rw') = —La{rw') > 
in (ri,r3) by Lemma 12.51 This implies that lwa{r) — rw'{r) < for r G (ri,r3) and a 
contradiction. Therefore 

lJLaWa{r) > rw' {r) for r G (0,ro), (2-18) 

and as a direct consequence, ■Wq > in (0, tq). 

It remains to check that Wa{rQ) > 0. Fix sq G (0,ro) and define /q := •5o^'('5o)/'"^a('5o) and 
'^o ■= (loWa — rw'y{so). We first observe that, since sq G (0,ro), (j2.18|) guarantees that 
Iq G (0, iJ,a). We next show that /cq > 0. Indeed, if fco < 0, we have loWa{r) — rw'{r) > in 
a left-neighborhood of sq as loWa{so) — soUj'{so) = 0. Defining 

si :=sup{r G (0,So) : loVUair) - ruj'{r) = 0} , 

we have just established that si < sq while Lemma 12.51 guarantees that si > as /q < M^^- 
Therefore, setting h{r) := loWa{r) —rw'{r) for r G [si,So], we have /i(si) = /i(so) = 
and La{loWa — rw') > in (si,so). This implies that h < in {si,sq) by Lemma F2.6I 
and a contradiction. Next, if /cq = 0, it follows from Lemma 12.51 that {p — l)sQ{loWa — 
rw')"{so) = La{loWa — rvu'){so) > from which we deduce that loWa{r) — rw'{r) > in a 
left-neighborhood of sq. Arguing as in the previous case, we are again led to a contradiction 
and conclude that /cq > 0. 

Consider now the solution ip of La{ip) = in (0,i?(a)) with initial conditions iI'^sq) = 
and tp'iso) = 1. We claim that -0 cannot vanish in (sqj'^o]- Indeed, since La{wa) = with 
Waiso) > sqvu' (sq) / fia > by Lemma [231 and ()2.18p . vua and ^p are linearly independent 
solutions to the same second-order ordinary differential equation and the oscillation the- 
orem guarantees that there is at least one zero of Wa between two zeros of ■0- Thus, if 
tpis) = ijj(so) = for some s G {sQ,ro], the function vUa has to vanish at least once in (sq, s) 
which contradicts (|2.18p . Therefore, tplr) > for r G (soj^o]- 

Define now the function ip := loWa — kQip. Then, by construction and Lemma l2.5[ 

La{ip - rw') > in (0, R{a)), ip{sQ) - sqw'{sq) = {f - rw')'{sQ) = 0, 

so that {p — 1)sq((/? — rw')"{sQ) = La(y' — fw'){sQ) > 0. Consequently, c^(r) — rw'{r) > in 
a right-neighborhood of sq. Assume for contradiction that there is 52 G (so,ro] such that 
Lp{r) — rw'{r) > for r G (sq, S2) and ^{32) — S2w'{s2) = 0. As La(v3 — rw') > in (sq, ^2), 
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we deduce from Lemma 12.61 that (p{r) — rw'{r) < for r G (so,S2) and a contradiction. 
Consequently, ip{r) — rw'{r) > for r € (sc^o]. In particular, 

loWairo) = 99(ro) + koip{ro) > row'{ro) + koTp{ro) > 0, 
which ends the proof. □ 
As in [6l [22], we split the range of a in three disjoint sets: 

A := {a > : there exists Ri{a) £ {0,R{a)), w'{Ri{a);a)=0}, 
B := < a > : w'{.; a) > in (0, oo), lim w{r; o) < oo > , 

C := < a > : w'{.] a) > in (0, oo), lim w(r; a) = oo\ . 

Since w'{r;a) ~ fiar^~^ as r — ?> 0, we have w'[.;a) > in a right-neighborhood of r = 
from which we deduce that the three sets are disjoint and Au BUC = (0, oo). According 
to the discussion at the beginning of Section 12.21 and the definition of w, a fast orbit 
corresponds to a € B while a slow orbit starts from a € C. Thus, the first step is to 
show that B is non-empty which will follow from the fact that A and C are both open 
intervals. In a second step, we shall prove that B reduces to a single point, thereby showing 
the existence and uniqueness of a radially symmetric self-similar very singular solution a 
stated in Theorem ll.il 

2.3 Characterization of A 

We show that, for all a £ A, the profile /(.;a) crosses the r-axis, that is, R{a) < oo. To 
this end, we introduce 



which is the constant solution of ()2.15p without the term r^\rw'{r) — fiw{r)\'^~^'^~^ which 
is expected to be negligible in the limit r — > oo as ?] < 0. 

Lemma 2.8. Let a > 0. Then, a £ A if and only if there exists Ri{a) € (0, R{a)) such that 
w'{.;a) > in (0,i?i(a)), w'{.;a) < in {Ri{a), R{a)) and w"{Ri{a)]a) < 0. Moreover, if 
a £ A, then 

sup w{r;a) < w* . (2.20) 

re{0,R(a)) 

Proof. We adapt the technique in [S]. Let a £ A and denote the smallest positive zero of w' 
in {0,R{a)) by Ri{a), its existence being guaranteed as a € A. Then, w' > in (0,i?i(a)) 
and w" {Ri{a)) < 0. Assume for contradiction that ■w"{Ri{a)) = 0. Differentiating (|2.15|) 
with respect to r and taking r = Ri{a) in the resulting identity give 

(p - l)Ri{afw"'{Ri{a)) = i]Ri{a)'^"^{fiw{Ri{a)f+'^~P < 0. 

Consequently, there exists some 6 > such that w" > in (i?i(a) — 6,Ri{a)), which 
entails that w'{r) < w'{Ri{a)) = for r G (i?i(a) — 5,Ri{a)) and a contradiction with the 
definition of i?i(a). Hence, w"{Ri{a)) < 0. 

We show next that ■w' < in {Ri{a) , R{a)) . Since w'{Ri{a)) = and w"{Ri{a)) < 0, we 
clearly have w' < in a right-neighborhood of Ri{a). Assume for contradiction that there 
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is i?2(a) G {Ri{a),R{a)) such that w' < in {Ri{a), R2{a)) and w'{R2{a)) = 0. Then 
w"{R2{a)) > 0. We evaluate ([213]) at Ri{a), i = 1,2, and find 

/x(^ - N)w{Ri{a)) + (a - /3/i)(/iu;(i?i(a)))2-Pzz;(i?,(a)) - Ri{ay {i^w{Ri{a))y-P+^ 
= -{p-l)Ri{afw"{Ri{a)), 

for z = 1,2. Since w"{Ri{a)) < < w"{R2{a)) we deduce 

ifiw{R,{a))f~P [/3m - a + /.i?i(a)''(/i^(i2i(a)))'^-i] < /x(a. - N) 

< {fiw{R2{a))f-P [/3/i - Q + f,R2{a)\i^w{R2{a))y~'] . 

But this is a contradiction, since 'w{Ri{a)) > w{R2{a)), 2—p > 0, 2—p+q—l = q—p+1 > 0, 
/3/i - a > and i?i(a)'' > i?2(a)''- Thus, w' < in {Ri{a), R{a)). 

The converse statement being obviously true, it remains to show ()2.20p . If a € A, w{.) 
has a single critical point Ri{a) in (0, R{a)) and u; attains its maximum at this point. The 
inequality ()2.20p then readily follows from the first inequality in ()2.2ip . □ 

The next step is to show that any profile /(.;a) starting from a & A crosses the r-axis. 
The following lemma which is implicitly contained in |21 Lemma 3.1] will be useful. 

Lemma 2.9. Let h be a nonnegative function in C^([0, oo)) such that there is a sequence 
{'^k)k>i, '"A; ^ oo as k ^ cx), for which h{rk) — > as k ^- cxd. Then, there is a sequence 
{Pk)k>i, Pfc — ^ oo as k ^ oo, such that h^p^) — > and Pkh'{pk) — > as k ^ oo. 

Proof. If h oscillates infinitely many times, then one can select {pk)k>i to be a sequence 
of local minimum points of h. Otherwise, if h does not oscillate infinitely many times, 
it eventually decreases monotonically to zero and h' belongs to L^(0,oo). Then, one can 
select pk > k such that pk \h'{pk)\ < 1/^- D 

The next result provides a complete description of A in terms of R{a). 
Lemma 2.10. Let a > 0. Then a (z A if and only if R{a) < oo. 

Proof. The proof is divided into two technical steps. Let a (^ A and suppose for contra- 
diction that R{a) = oo. 

Step 1. Interpretation in terms of decay. By Lemma 12.81 we deduce that w' < in 
(i?i(a), oo) and the non- negativity of w guarantees that w has a limit ^ > as r ^ oo. In 
addition, w' belongs to L^(i?i(a), oo) from which we deduce that there exists a sequence 
Tfc — >■ oo such that rkw'{rk) — > as A; ^> oo. Since w' < 0, it follows from Lemma [2.91 that 
we may also assume that rk{rkw"{rk) + w'{rk)) — > as /c — )■ oo, whence r1w"{rk) — > 
as /c — )■ oo. We evaluate (j2.15p at r = r^ and pass to the limit as fc ^ oo to deduce that 
i € {0, w*} (recall that 77 < 0). Owing to ()2.20p . the possibility i = w* is excluded and we 
conclude that 

lim w{r) = 0. (2.22) 

Step 2. The contradiction. Since p < 2, we infer from (|2.14p and ()2.22p that 
^-^^^ < C rf{r)^^-PyP = C u;(r)(2-P)/P -^ 0. 
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Consequently, there exists r* > Ri (a) such that 

rf'ir) 
— u < „, , < tor any r > r^,, 
fir) 

from which we deduce that 

w'{r) = ri" {f'{r) + ^xf{r)) > , r G (r.,, oo) , 

and a contradiction with (j2.22p . Therefore R{a) < oo. 

Conversely, if R{a) < cx), then w{R{a)) = = w{0), which implies that w has a maximum 
point in (0, i?(a)), hence a G A. Q 

Proposition 2.11. The set A is an open interval of the form (0,a*) for some a* > 0. 

Proof. We first show that A is non-empty and that (0, a) C A for a sufficiently small. The 
proof is the same as that of [221 Theorem 2] and we only sketch it here for the sake of 
completeness. We perform a rescaling in the r-variable and, for a > 0, we define g{.; a) by 

/(r; a) = ag ('ra^^-rf/P; a\ , r € [0, R{a)) . 
Then g = g{.;a) solves 



iWr^gJis) + ^^^(|5r"V)(s) + I3sg'{s) + ag{s) - a(2'?-rt/P|^'(,)|</ = q. 



(2.23) 



[ 5(0) = 1, 5'(0) = 0. 
Owing to Lemma |2.H we have 

5(s)>0 and - a^/« a-^^g-rt/pg < ^'(g) < q for s G fo, i?(a)a(2-P)/p) . (2.24) 



We next study the limit of (|2.23p as a ^- which reads, since q > p/2, 

(|/i'|P-2/i')'(s) + ^^^^{\h'\P~^h'){s) + /3s/i'(s) + ah{s) = 0, 
s 

h{0) = 1, h'{0) = 0. 



(2.25) 



Arguing as in Lemma 12.11 there is Sq > such that the solution h to (|2.25p satisfies 
/i(S'o) = 0, h{s) > and h'{s) < for s G (0, Sq). Moreover, proceeding as in the proof of 
[22t Theorem 2], we show that 5o < oo with h'{So) < 0. It is then easy to show, using the 
continuous dependence of solutions for (|2.23p . that for a > sufficiently small, g{.; a) also 
vanishes and thus R{a) < oo. Consequently, thanks to Lemma [2. 101 (0, a) C A for a small 
enough. 

It remains to show that A is an open interval. By Lemma 12. 8| if a € ^, we have 
w'{Ri{a);a) = with w"{Ri{a);a) < 0. By continuous dependence, this property read- 
ily implies that A is open. In addition, we can apply the implicit function theorem and 
conclude that a i->- Ri{a) belongs to C^{A). Define m{a) := w{Ri{a);a) for a G A. Then 
m G C^ (A) and we infer from Lemma 12.71 (with tq = Ri (a) ) and Lemma 12.81 that 

-^(a) = w'(Ri(a);a)—^(a) + daw(Ri(a);a) = daw(Ri(a);a)>0, aeA. (2.26) 
da da 
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Consider now a2 (z A and define ai > by 

oi := sup{a e (0,02) : a ^ A} . 
Since A is open, we have ai < 02, ai A, and (01,02) C ^. Setting 

p := lim infill (a) G [0, 00] , 

there are three possibihties: 

If p = 00, we actuahy have Ri{a) — > 00 as a \ ai and it follows by continuous de- 
pendence that, for any r > 0, we have r < Ri{a) for a > ai close enough to ai and 
thus 

w'{r;ai) = lim w'{r;a) > 0. 

a\ai 

In addition, by Lemma 12.101 and ()2.26p . 

w{r; ai) = lim w(r; a) < ?7i((ai + 0,2)/'^) < w*, (2.27) 

a\ai 

SO that w{.; oi) is a non-decreasing and bounded function. It thus has a limit i G [0, m((ai + 
02)72)] as r — )• 00. If £ > 0, we argue as in Step 1 of the proof of Lemma 12.101 to obtain 
that £ = w* and a contradiction. Consequently, £ = so that w{.; ai) = and thus ai = 0. 

If /? € (0, 00), there is a sequence (a^) € (ai, 02) such that a^ — > ai and Ri{a^) — > p as 
j -^ 00. By continuous dependence we have 

w'{p;ai)= lira w'{Ri{a-');a-') = , 

j-s-oo 

whence ai (^ A and a contradiction. 

If p = 0, there is a sequence (a-') G (01,02) such that a^ — > oi and Ri{a^) — > as 
j -^ 00. Let us assume for contradiction that oi > 0. Then there is ri > such that 
t(;'(r; oi) > for r G (0, ri) by Lemma [2^21 Given r G (0, ri), we have Ri{a^) < n for j 
large enough whence 

w'{r;ai) = lim w'{r;a^) < 0, 

by continuous dependence, and a contradiction. Consequently, oi = in this case as well. 
We have thus shown that, given 02 G A, the interval (0,02) is included in A from which 
Proposition 12.111 follows. [] 

2.4 Characterization of C 

We begin with the following useful result. 

Lemma 2.12. Let o > 0. Then a ^ C if and only if 

sup w{r]a)>w*. (2.28) 

re(0,iJ{a)) 

Proof. The direct implication is obvious. Conversely, if o > is such that ()2.28p is satisfied, 
it follows from Lemma [2.8l that a ^ A. Then o G BUC and w{.; a) is an increasing function 
in (0, 00). If it is bounded, it has a finite limit ^ as r ^- 00 and we argue as in the first 
step of the proof of Lemma [2.101 to deduce that i G {0, u;*}, clearly contradicting ()2.28p . 
Consequently, w{.;a) is unbounded and a G C. Q 
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Proposition 2.13. The set C is an open interval of the form (a*,oo). 

Proof. Let us first show that C is non-empty. Given a > 0, by Lemma 12.11 we have 
/'(r) > -(aa)^/"? for r G {0,R{a)), whence 

/(r) >a- {aa)^/''r for r E (0, /?(a)). (2.29) 

In particular, R{a) > 0(00)-^/" = a-'^^'^a'^i-^y" . Consequently, a-'^/ia'-'^'^yi/l G (0,i?(a)) 
and it follows from (|2.29p that 



w 



(a-i/'^a('?-i)/V2;a) > -(«-i/'^a(''"i)/V2)^ = a^^/^^^-p) 2-2/(2-^) ^[2g~p) M2-p) ^ 



which can be chosen greater that 2w* for a large enough since q > p/2, and guarantees that 
a (z C for a large enough, according to Lemma l2.12[ This shows that C is non-empty and 
that a € C for a large enough. In addition, C is clearly an open set thanks to Lemma [2. 121 
Now, let ao G C. Then w'{.;ao) > in (0, cxj) and also daw{.;ao) > in (0,oo) by 
Lemma 12.71 Therefore, C is an open interval. □ 

We now show that the set C is coinposed in fact of the slow orbits. 
Proposition 2.14. For every a (z C, there exists k{a) S (0,oo) such that 

lim AP-''yi'>-P+^')f{r;a) = k{a). 

r— >oo 

While the outcome of Proposition 12.141 is similar to that of [6l Theorem 4.1], only the 
first step of the proof of Proposition 12.141 borrows some ideas of the proof of [6l Theorem 
4.1], that is, the introduction of a new variable and a new unknown function. The other 
two steps are different owing to the different nature of the absorption on the one hand 
(Step 2) and the fact that q takes values below 1 on the other hand (Step 3). Concerning 
the latter, we prove a general lemma. 

Lemma 2.15. Let X € C^([0, 00)) satisfy the differential inequality 

X'(i) + C5eT*X(t) <C76(l + e'5*), i>0, (2.30) 

for some < 5 < 7 and C5 > 0, Cq > 0. Given 9 G (0,7 — 6), there exist tg > and 
Ke > depending only on 7, S, C5, Cq, ^(0), and 9 such that X{t) < Kee~^^ for all 
t>te. 

Proof. Make the change of variable X{t) = Y{e'^^) = Y{t) for t > 0, where r := e"^* > 1. 
Then the inequality (|2.30p becomes 

7y'(r) + C5y(r) < ^ (l + r-^/^), r > 1, 
or, equivalently, after a straightforward transformation 

A(KMe--)<Co(i + ^)e-". r>l. 

with C7 := C5/7. After integration over (1,t), r > 1, we obtain 

Y{t) < y(l)e-^^(--^) + C, f^ (i + ^(:^) e^^('^--) da, r > 1. (2.31) 
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We take r > 2 sufficiently large such that r/ log r > 2/Cj. In order to estimate the integral 
in the right-hand side of (j2.3ip . we split it and handle differently the contributions to the 
integral of the intervals (l,r/2), (r/2,r — (logT)/C7) and (r — (logr)/C7, r). First, 

Next, a -T < -(logT)/C7 for a € (r/2,r - (logr)/C7) so that 

-{logT)/C7 /I \ \ ^, , 1 /■^-{log^)/C7 /I 1 \ 

- + ^-^T^ e^^("-") da < - - + , ,,, da 



cr o-(t~^)/t/ ~ T L/9 IcT cr(T-'^)/')' 



t/2 V" O-^' ^' '/ ' Jt/2 

o-=t/2 



1 r 
< - 

T 



log (J + - cr 




!„■•/']'" 



< l('-+ ' 



5 yr r(T~^)/T, 
Finally, if o" S (r — (log t)/C7, r), we have a > t/2 and a — r < and 

7r-(logr)/C7 Vf^ C7(^-'5)/7y yr-(logr)/C7 V^ V^/ / 



do" 



^ 21ogT /l^ 1 



C7 \t T(i-&)h 

Plugging these estimates into (j2.3ip . we obtain that, given q € (0, (7 — 5)/5), there are 
Kg > and Tg > 2 such that Y{t) < K^t^^ for all r > r^. Coming back to X gives the 
claim. □ 

The proof of Proposition 12.141 is divided into several steps. 

Proof of Proposition \2. 14\ Step 1. Introducing a new function. Fix a £ C. As in [6l 

Theorem 4.1], we define a new function A by 

/(e^;a) = /(l;a)exp('-y""A(s)dsV r e M. (2.32) 

Then 

/'(e-;a) = -e-^A(T)/(e^;a), /"(e^;a) = e'^^ f {e^ ; a)[A{T) + Airf - A'(r)], 
and it follows from (I2.15P that A solves the following differential equation 

(p-l)A'(T) = F(r,A(T)), TGM, (2.33) 

with 

F{t, A) := {p - 1)A^ + {p-N) A + A^~Pw{e^;af~P [a - /3A - A^'e''^^ f{e^;ay~^] . 

Since a G C, we have R{a) = 00 and it follows from Lemma l2. II that f'{.;a) < 0, whence 
A > 0. In addition, for r € M, 

w'{e^; a) = e^^"'^- (eV'(e"; a) + /x/(e-; a)) = e(^"i)-(;U - A)/(e-; a) > 
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since a ^ C, from which we deduce that A(t) < /i for all r G M. 
Step 2. Limit of A as r ^ oo. We claim that 

lim e-''^/(e^;a)'?-^ =0. (2.34) 

T— >00 

Indeed, either q > I and (|2.34|) is obvious thanks to the boundedness of / which guarantees 
that 

or g S (p/2, 1) and 

and both terms converge to zero as r — t- oo since a & C and q € (p/2, 1). Hence ()2.34p 
holds true. Consequently, since a G C, we have 

a OL 

lim F(r, 7) = +cx3 if 7 < — , lim F(r, 7) = —00 if 7 > — , 

T— 7>oo p r— >oo p 

which ensures that 

lim A(r) = -. (2.35) 

T^OO p 

Step 3. Exponential decay of A — a//5. We study the rate of convergence of A(r) to 
a//3 as r — > 00. This last step of the proof is done in [6] by construction of suitable sub- 
and supersolutions. We give here an alternative and direct approach based on Lemma [2. 151 
which allows us to handle also the case q <1. 
Since A(r) — > a/ (3 as r — ?• 00, taking into account the definition of A, we find 

rf'ir) a 

hm ^-Y = --, 

hence, for any e G (0, a/ (3)^ there exists r^ > 1 such that 

- ^ -e < ^7^ = -A(logr) <--+£, r>re. (2.36) 

By integrating, we obtain that there exists a constant Kf, > such that 

K^r-^^/P)-e < j(^) < /^^r■-("/^)+^ r > r,. (2.37) 

We then come back to the differential equation (j2.33p and check that the conditions re- 
quired to apply Lemma 12.151 with either X = (A — a//3)+ or X = (a//5 — A)^ are fulfilled. 
Indeed, we can write (j2.33p in the form 

{p - l)A'(r) + l3A{Tf-Pwie-;af-P (^A{t) - ^) = 5(r) , r € R , (2.38) 

with 

S{t) := {p - l)A(r)2 + (p - iV)A(r) - A{Tf-P+^w{e^; af-Pe-'^^ /{e^; af-^ . 

On the one hand, owing to the boundedness of A, the non-negativity of /(.; a), w{.; a), and 
A, and (|2.37p . we have 

S{t)<{p-1)h^+PH, tGM, (2.39) 
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and 

S{t) > -Nfi - /i«-P+2 e^P-"^^ fie-"; 0)"-^+^ 

> —Nil - uF~'P^'^ gCp-?)^ ^g-p+i g(e(9-p+i)+'?-p)T 
S{t) > -Nfx - /i«-P+2 K9-P+1 e^(g-P+i)- , r > logr, . (2.40) 

On the other hand, it follows from ([235]) . ([236]) . and (fOT]) that 

M(r)^-Pu;(e^;a)2-P>/3 (^-e) ^2-^ g(2-p){/.-{a//3)-e)r ^ T > log r, . (2.41) 

We now fix e € (0, a//3) such that 

a 1 (01 

^-28 ^""^ 2e(g-p+l) < --e(2-p) = (2-p) U---e 

such a choice being always possible since /i > a//?. Introducing Xi := (A — a//3)_|_, we 
infer from (12381) and (1239)1 that 

(p- l)X;(r) +/3A(r)2-Pzz;(e-;a)2-P Xi(r) < 5(r)+ < (p- 1)/^^ +p^, 

while (|2.4ip . the non- negativity of Xi, and the choice of e guarantee that 



{p-l)X[{T) + [i (^y ' e(2-P)(^-("/«-)^Xi(r)<5(r)+<(p-l)/.2+; 



for T > logr^. Consequently, Xi satisfies a differential inequality of the form of ()2.30p 
with 7 = (^ — {a/ 13) — e){2 — p) > and 5 = 0. We then apply Lemma 12.151 and conclude 
that 



(a{t) -^^ <C e-^^^-P+i)^ 



for T large enough. Similarly, setting X2 ■= (a//3 — A)^-, it follows from ()2.38p and ()2.40p 
that 

(p-l)X^(r) + /3A(r)2-Pu;(e-;a)2-PX2(r) < (-5(r))+ 

and from (|2.4ip , the non- negativity of Xi , and the choice of e that 

(p-l)X^(r)+/5 (^)'~' e(2-^)('^-("/^)-)^ X2(r) 

for T > logr^. Consequently, X2 satisfies a differential inequality of the form of (j2.30p 
with 7 = (/x — {a/ 13) — e)(2 — p) and (5 = e((7 — p + 1), and the choice of e guarantees that 
7 — (5 > (5 > 0. We then apply Lemma 12.151 and conclude that 



(^ - K{t)\ < C e-^(9-P+i)- 
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for r large enough. We have thus estabhshed that |A(r) — q//3| < C e ^'"^ p+i)"^ for r large 
enough. In particular, A — a/ (3 belongs to L^{1, oo) and, recalling ()2.32p . we realize that 

r"/^/(r; a) = /(I; a) exp (- £°' ' (A(r) - ^) dr^ 

_^/(l;a)exp (^-^" (^A(r) - ^^ dr) =: A:(a) G (0, oo) , 
as stated in Proposition 12.141 □ 

2.5 Characterization of B. Uniqueness 

Since A = (0, a*) and C = (a*,oo), we already know that B is non-empty, that is, there 
exists at least one radially symmetric self-similar very singular solution. We complete the 
proof of Theorem 11.11 bv the uniqueness result. 

Proposition 2.16. We have a* = a* , thus B contains only one element. 

Proof. Let a £ B. Then, R{a) = oo, w{.; a) is increasing and has a finite limit as r ^ oo. 
Arguing as in the first step of the proof of Lemma I2.10| we deduce that 

lim r^f{r\a) = lim w{r;a) = w* . (2.42) 

r— >oo ■; — ^oo 

As in [22], the proof is divided into two steps, with a different argument to prove the first 

step. 

Step 1. We show that rw'{r;a) — > as r ^ oo for a £ B. Indeed, we notice that the 

differential equation ()2.ip also reads 



|:(r^- V'Wr V'(r) + /3r^/(r)) = r^-^l/'Wr " (« " m)f{r)) , r > . 
Since q > p/2, it follows from Lemma 1 2. II and Proposition 12.41 that, as r — t- oo, 

|/'(r;a)|''<C4/(r;a)2''/P = o(/(r;a)), 
and we deduce from (|2.42p that, as r — > oo. 

Since fj, > N, we obtain after integration 



-r^-\-f\r;a)f-^ + pr^ fir;a) r. [a - l3N)w*- 



hence, using again (j2.42p . 

-i-fXr; a)r-' ~ -/3r/(r; a) + ^^^^-'^^'^^ - (^^ " /?) -^-'^^"^^ 



as r ^- oo. Noticing that 

a- f3N 



P = - ... , .. ... <o, 



fi-N "^ p{N + l)-2N 
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we conclude that, as r ^^ oo, 



fir; a) 



w 



i/(p-i) 



-2/(2-p) 



^p{N + l)-2N^ 
Consequently, 

lim rw'{r;a)= lim (r^/^^~P^ f'{r;a) + firP/^^~P^ f{r;a)) 



^ 1/(P-1) 






^p{N + l)-2N 
But from the definition and a simple computation we find 

w* = ^i^P'^'^l^'^-P\p{N + 1) - 2iV)i/(2-p)^ 

that is 

^ = (w;*)(2-p)/(p-i)(p(Ar + 1) - 2Ny^l^P~^\ 

hence 

i/(p-i) 



/.u;* = (u;*)i+(2-P)/(P-i)(p(Ar + 1) _ 2iV)"i/(P-i) = l- 

\1 



w 



^p{N + l)-2N^ 

We have thus shown that 

lim rw'{r;a) =0. (2.43) 

r— >oo 

Step 2. We now follow the argument in [51 Theorem 5.1]. Let a ^ B. Since w{r; a) — > w* 
and rw'{r;a) — > as r ^ oo by ()2.42p and ()2.43p . the homogeneous linear operator La 
defined by (|2.16|) can be written as 

La{^) = r\"{r) + Pi(r) r (^'(r) + Po{r) (^(r) , 

with 

, . , ^-1 p{N + 1) -2N , 

lim Pi r = h := - 2/i + ^^ ^ \ 1 + /3/i , 

P-I P-I .2 44) 

p{N + l)-2N ^ ' 

lim Po{r) = -c , c := ^— ^- > . 

r— >oo r) — 1 

Owing to the negativity of c, we claim that, given two linearly independent solutions of 
the equation La{(p) = 0, one of them has to be unbounded as r — )• oo. Taking this result 
for granted, we define the function ip as the solution to La{ip) = in (1, oo) with ^(1) = 0, 
'(/''(I) = 1. Arguing as in Lemma 12.71 with rg = oo (recall that R{a) = oo as a € B), we 
conclude that there is a positive real number k > such that daw{.; a) > Ktp > in (1, oo). 
Since La{daw{.; a)) = with daw{l;a) > by Lemma [231 and Lemma \T7\ daw{.;a) and 
■0 are clearly linearly independent solutions of the equation La{f) = 0. The above claim 
implies that one of them must be unbounded and thus that daw{.;a) is unbounded since 
it dominates ip. 

We then deduce that, for all a G B, daw{r; a) — > oo as r ^^ oo. Therefore, if a* < a*, we 
can apply Fatou's lemma to get 

a* a* 

= lim [w{r,a*) — w{r,a^,)] = lim / daw(r; a) da > / lim. inf daW^r; a) da = oo, 

r— >oo r— >oo J J r— >oo 
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which is a contradiction. Hence a* = a* . 

Let us finally sketch the proof of the above claim. Introducing the new variable r = log r, 
the linear operator La transforms into the linear operator 

LaWir) := (^"(r) + A(r) ^'{r) + Po{t) <^(t) 

with Pi{t) := Pi(e^) — 1 and Pq{t) := Po(e^)- Introducing next the change of function 

V'(t) = exp Q y Pi(s) ds] (^(r) , r > , 

we obtain the canonical form Ka of the operator La which reads 

52 p' 

KaWir) := r{r) + Q{t) i;{T) with Q -.= P^ - ^ ^. 

At this point, we note that (j2.44p ensures that 

52 
lim Q{t) = -A§ := -c - — < . (2.45) 

T— >-00 4 

Then, if ip is a solution to Ka{ip) = 0, the variation of constants formula implies that there 
are V'o ^ I^ and -01 G M such that 

V,(r) = V'o e-^o- + V'l e^o- + T j^iiMfl (e^o(.-.) _ g-Aol.-.)) | rf. 



,^ ^ ,. „., r>0. 

I ^^0 



Arguing by contradiction, it is now easy to show that, if ^i > and A G (0, Aq), then the 
function r 1-^ e~^'^'4){T) cannot be bounded. Coming back to La gives the expected result 
after choosing A appropriately close to Aq. D 



A Viscosity solutions to (11.11 ) 



We recall here the definition of viscosity solutions for singular diffusion equations such as 
(jl.ip as introduced in |lltll7|. Owing to the singularity of the diffusion, it differs from the 
usual definition of vicosity solutions (see, e.g., [3) by restricting the class of comparison 
functions. More precisely, let H be the set of functions ^ € C2([0, oo)) satisfying 

^(0) = e'(0) = ^"(0) = 0, i"ir) > for ah r > 0, lim W {rW-"^ i" (r) = 0. (A.l) 

■; — s-0 

For example, ^{r) = r" with a > p/{p — 1) > 2 belongs to H. In fact, if ^ € H, it readily 
follows from (|A.ip that 

^(r) 

We next introduce the class A of admissible comparison functions ip S C^iQ oo) defined 
as follows: ^ € .4, if, for any {tQ,XQ) S Q^o where Vip{tQ,xo) = 0, there exist a constant 
(5 > 0, a function ^ £ H, and a nonnegative function u £ C([0, oo)) satisfying uj{r)/r — > 
as r ^- 0, such that, for all (t, x) G Qoo with jx — xo| + |t — to| < 5, we have 

|V(t, x) - ^P{to,xo) - dttP{to,xo){t -to)\< C{\x -xo\)+ uj{\t - tol). (A.3) 

21 



Definition A.l. An upper semicontinuous function u : Qoo — > M is a viscosity subsolution 
to (jl.ip in Qoo if, whenever ip € A and {tQ,xo) G Qoo o.'re such that 

u{to, xo) = ipito, xo), u{t, x) < ipit, x), for all {t, x) G Qoo \ {{to,xo)}, 

then 

dtipito,xo) < Ap^{to,xo) -\V^p{to,xo)\'^ ifVip{to,xo)y^O, .. . 

dMto,XQ)<0 i/VV(to,xo) = 0. ^ ■ ' 

A lower semicontinuous function u : Qoo -^^ is a viscosity supersolution to (jl.ip in Qoo 
if —u is a viscosity subsolution to (jl.ip in Qoo- A continuous function u : Qoo — )> M is a 
viscosity solution to (jl.ip in Qoo if it is a viscosity subsolution and supersolution. 
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